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formulae. An accuracy gain can be achieved even for the situation when two of the four 
orbits are ghosts. This is illustrated for a model system, the kicked top. 
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1. Introduction 

Periodic orbits provide the skeleton of the dynamics of classical Hamiltonian systems. 
Generic dynamical systems display a mixed phase space, consisting of islands of stability 
residing in chaotic seas. The periodic orbits are neither grouped in families, like in integrable 
systems, nor are they well isolated and all unstable, as for chaotic systems. One rather finds 
also stable periodic orbits surrounded by islands of regular behaviour. These islands of 
stability look locally like an almost integrable system, altogether with a KAM structure 
of invariant tori and chains of periodic orbits, remnants of rational tori of a supposedly 
contiguous integrable situation 

In a semiclassical treatment of the corresponding quantum system, clusters of proximate 
orbits display a collective behaviour. The bifurcations at the centre of the island and their 
semiclassical treatment have been addressed in a number of recent works, both for the generic 
variants p|, ^, ^, ^ ^ as well as for classically non-generic, but semiclassically still relevant 
cases ^. 

In H| the class of near-integrable systems has been addressed, and a uniform semiclassical 
approximation for the most frequently encountered broken rational tori (consisting of a stable 
and an unstable periodic orbit) was presented. We will call these tori the 'simple' tori. 
Interestingly, the semiclassical approximation works reasonably well even beyond the point 



where the stable orbit becomes unstable |10 . The same configuration of a stable and an 



unstable orbit is also typical close (in parameter space) to most types of period-n-tupling 
bifurcations at the centre of a stability island. For bifurcation number n > 5, for the island- 
chain scenario with n = 4, and also (as a consequence of a more complicated bifurcation 
scenario) for n = 3, two satellite orbits are expelled from the centre. At a certain distance 
from the bifurcation the satellite orbits can be treated as isolated from the central orbit; 
however, a collective semiclassical treatment of the two satellites is often still necessary, and 
this can be achieved by using the abovementioned approximation for the simple torus. 

Although encountered less frequently, there are situations where a broken torus consists 
not of two, but a higher number of periodic orbits. In the islands of stability tori of this 
type appear especially at larger distance from a bifurcation. In this work we study the 
twofold-broken rational torus, consisting of two stable and two unstable periodic orbits. It 
is described by a normal form which is obtained from the normal form of the simple torus 
by including the second harmonic in an angular coordinate. From the normal form we 
construct a uniform approximation that can be used to improve semiclassical trace formulae. 
Indeed, the relevance of this configuration is much enhanced in the semiclassical context: 
Here one has to consider also 'ghost' orbits with complex coordinates and even a simple 
broken torus can be affected by nearby ghosts, making a treatment as a 'pre-formed' twofold- 
broken torus advisable. This is illustrated in a model system, a periodically driven angular 
momentum vector (the kicked top), where we find a configuration of four period-three orbits 
which can be regarded as a twofold-broken rational torus. A reduction of the error of the 
trace formula by a factor of about 2 — 3 is found even when two of these satellites are ghosts. 
An even higher accuracy gain is attained for so-called 'inverse- fi, spectroscopy'. 



2. Normal forms and uniform approximations 

We restrict the analysis to two-dimensional area-preserving maps. (The results are also 
applicable to autonomous Hamiltonian systems with two degrees of freedom.) The quantum 
version of the map is generated by the unitary Floquet operator F which acts on the vectors 
of a Hilbert space, mapping the space onto itself. The semiclassical trace formula relates the 
traces trF" to the classical periodic orbits. Isolated orbits of primitive period no give an 
additive contribution ||, H |l|, |l| 



C = A exp 



(1) 
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with amplitude 

A = no|2-trAf|"i/2 



(2) 



to all traces tr with n = uqt and integer repetition number r. Besides the primitive 
period, three classical quantities of the (r-th return of the) periodic orbit enter, the action 
S, the trace of the linearized n-step map M , and the Maslov index /i. 

The expression (|l|) is derived by a stationary-phase approximation and becomes 
inaccurate when orbits lie close together. Then a collective treatment of the region 
inhabited by the proximate orbits becomes necessary. This can be achieved by introducing 
normal forms for a phase function $ and an amplitude function 5" into the more general 
expression ||, Q 



Cn= 2i^d^' d/vl/(y,',/)exp 



(3) 



Here /, ip are canonical polar (or cylinder) coordinates, and v is the Morse index. 

The famous Poincare-Birkhoff theorem states that a perturbation of an integrable system 
causes tori with rational winding number to break into chains of alternating stable and 
unstable periodic points. However, it does not give a quantitative criterion for the number 
of distinct orbits that lie on this chain. The simple broken torus is described by the normal 
form 



S{I, ip') = So + Itp' -al^ -b cos if' 

for the generating function S, with constants Sq, 
(/, ifi) — > (/', (f'), implicitly given by 



dS_ 



I' = 



dS 
97' 



(4) 

and h. The corresponding map 



(5) 



is the well-known standard map. This normal form has been used in Q to obtain a uniform 
semiclassical approximation for the simple broken torus, smoothly interpolating between the 
two non-commuting classical and integrable limits (?i — » and & — > 0, respectively). 

A more complete picture can be obtained when one includes the second harmonic in the 
angular variable and works with the extended normal form 



5(1, </?') = Sq + lip' - al'^ - 6cos</?' - ccos{2ip' + 2ipo)- 



(6) 



The strategy of including higher-order terms in normal forms has been pursued before, 
with two different incentives. Firstly, the inclusion of higher orders is a tool to equip a 
normal form with a sufficient number of independent parameters. This can be necessary to 
account for all classical properties (stabilities and actions) of the periodic orbits described 
by the normal form. Secondly, the higher orders describe additional periodic orbits, and 
hence more complicated configurations than the usual normal forms |0, ||, |l^, A 
semiclassical description often succeeds only when all orbits of an extended normal form are 
treated collectively. Presently we aim at the inclusion of additional periodic orbits. 
The periodic orbits satisfy the fixed point conditions 



/ — = / + bsini^' + 2csin(2(^' + 2(po) 



dS 
dip' 
dS 



2al, 



resulting in / = and 



6sin(^' + 2csin(2(y5' + 2(po) = 0. 



(7) 
(8) 

(9) 



This condition amounts to finding the roots of a fourth-order polynomial in sin Lp' . Depending 
on the parameters 6, c, and lpq there are either four real solutions or two real and two complex 
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solutions. For |6| < 2|c| there are always four real solutions. For |6| > |c| only two real 
solutions are found. When \b/c\ is fixed in the range (1, 2) and (f^ is varied one finds tangent 
bifurcations with two real solutions on one side of the bifurcation and four real solutions 
on the other side. The real solutions correspond to conventional periodic orbits while the 
complex solutions are 'ghosts'. They are of no consequence for the classical dynamics but 
can be important for the semiclassical description of the quantum system, as has been shown 
in jl^ and as we shall see once more below. 

The joint contribution of the orbits on the twofold-broken torus is found by introducing 
into Eq. (^) the normal form 

$ = S'o - a/2 - & cos ip' -c cos(2(^' + 2(Pq) (10) 

[cf. Eq. (H)] for the phase function and 

* = 1 + dcos{f' + ipi) + ecos(2(^' + 2ipo) (11) 

for the amplitude function. The stationary-phase limit of Eq. (^ is a sum of four additive 
contributions of form (^, each representing one orbit. The four parameters 5*0, 6, c, ipo are 
determined by matching the phases of each contribution to the actions S of the periodic 
orbits. The parameters a, d, e, and tpi are fixed by the stability amplitudes A. This strategy 
works also when two of the orbits are ghosts: Phases and amplitudes become then complex, 
but are related by complex conjugation, and the number of real independent parameters 
remains unchanged. 

Eqs. (|3|, [10| , Pt|) represent a uniform approximation of the joint contribution of the twofold- 
broken torus. The integration over / is readily carried out, which leaves us with a one- 
dimensional strongly oscillating integral over the coordinate ip' . Numerically it is most 
conveniently evaluated by the method of steepest descent, for which the integration contour is 
deformed into the complex plain. On the new contour the integrand decreases exponentially. 
The new contour can also visit stationary points with complex coordinates, i. e., ghost orbits. 



3. Numerical results 



We shall illustrate our findings, and especially the relevance of ghosts, for a situation 
encountered in a model system, the dynamics of a periodically driven angular momentum 
vector J (the kicked top |l^). The components of J obey the usual commutation rules 
[J'xiJy] — iJz (and cyclic permutations). The total angular momentum = j[j -I- 1) 
is conserved, restricting the dynamics to the irreducible representations of the angular- 
momentum algebra. The Hilbert space dimension is 2j 
is l/(j + 1/2) and the classical limit is attained for j - 
operator of the explicit form 



1. The effective Planck's constant 
oo. We work here with a Floquet 



F 



exp 



J' 



2j 



X exp 



i/c.T 



- 1 



iPzJz 



exp [-ipyJy] 



2j + l 



(12) 



The dynamics consists of a sequence of linear rotations by angles pi alternating with torsions 
of strength ki. We hold the pi fixed {px — 0.3, Py — 1.0, Pz — 0.8) while varying the control 
parameter k = — lOk^- Complete semiclassical spectra of this system throughout the 
full transition from integrable {k — 0) to well-developed chaotic behaviour (fc « 10) have 
been presented in |2^. The system has also been used to illustrate uniform semiclassical 
approximations for various kinds of bifurcations ^, ^ . 

We concentrate on a particular configuration of period-three orbits which suggests a 
treatment as a twofold-broken torus. The configuration comes about in a sequence of three 
bifurcations: At A: = 1.9715 ... a pair of period-three satellites is born in close vicinity to a 
period-one orbit in the centre of a stability island. At A; = 1.9753 a period-tripling bifurcation 
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Figure 1. Phase space portrait of the kicked top with k = 4.0. The phase space is 
the unit sphere, parameterized here with the azimuthal angle q and the ^-component p (in 
conventional spherical coordinates, q = (p and p = cos 6). The circle indicates the central 
period-one orbit, the other symbols indicate the positions of the four period-three satellites 
(all orbits are unstable). 





Re 9 Re 9 

Figure 2. Integration contours in the complex (/j'-plane on both sides of the tangent 
bifurcation in the kicked top at = 3.7856. For each k the normal-form coefficients of the 
phase function <& have been determined from the actions of the periodic orbits. The contours 
pass through the saddle points of the normal form (circles) and proceed in the direction of 
steepest descent (constant phase of <3>). 



takes place where the stable satellite collides with the central orbit. On the other side of 
the bifurcation the satellites form a simple broken torus around the centre. This sequence of 
bifurcations can be described by extended normal forms |j, and the close neighbourhood 
of the two bifurcations is not exceptional (see e.g. |2^). As the control parameter k is 
increased further the satellites move away from the centre. At fc = 3.7856 another pair of 
period-three satellites appears in a tangent bifurcation. For k < 5 the four satellites form a 
configuration that resembles a twofold-broken torus. Figure]^ displays a phase space portrait 
for fc = 4, and Figure |^ shows steepest-descent contours in the complex (^'-plane on both 
sides of the final bifurcation. 

We evaluated ti at fc = 3.0 with the quantum number j ranging from 1 to 50. At 
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Figure 3. Semiclassical error of tr for the kicked top (k = 3.0) as a function of j in the 
two approximations i), ii) which are explained in the text. 




Figure 4. Collective peak of the four satellites and the central orbit in the action spectrum 
|r(S')p for k = 3.0. The semiclassical approximations i) and ii) (see text) are compared to 
the result of an exact quantum-mechanical computation. 



the given value of the control parameter only one of the pairs of satellites mentioned above 
has real coordinates, and its distance to the centre of the stability island is already quite 
large. The other two satellites are still ghosts, but their bifurcation is not far away. This 
leaves us with the choice between two semiclassical approximations: i) We can group the two 
real satellites together with the central orbit and treat the ghosts separately, or ii) we can 
group the four satellites as a twofold-broken torus and treat the central orbit separately. (In 
principle, we could complicate matters even more and group all the orbits together, but this 
is rather impractical and beyond the scope of the present work.) The semiclassical evaluation 
of the trace involves also five other orbits. The error of the two approximations is shown 
in Figure ^. The error of approximation ii) is about a factor 2 — 3 smaller than that of 
approximation i). Although the accuracy gain is not dramatic, this result favours clearly a 
treatment of the satellites as a pre- formed twofold-broken torus. 

A somewhat more demanding test is aided by 'inverse-?! spectroscopy' [|[ |2|. We consider 
the discrete truncated Fourier analysis 

32 

^(^) = ^EtrF3(j)exphi(j + 1)5] (13) 
of the trace with respect to the quantum number j. The 'action spectrum' |r(S')p displays 
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peaks at the actions of the periodic orbits that contribute to trF^. Since accidental 
action degeneracies do not occur in the present example, the quality of the semiclassical 
approximations i) and ii) can now be assessed directly, without interference of the remaining 
orbits. Figure ^ shows the collective peak of the four satellites and the central orbit. 
Approximation ii) agrees almost perfectly with the exact result, while approximation i) 
overestimates the peak-height distinctively. 

Of course, the uniform approximation presented here is not restricted to the situation 
where two of the orbits are ghosts, but is valid for four real orbits on the broken torus as 
well. Indeed we find an improvement in the semiclassical accuracy of tr over the full range 
2.8 < fc < 5. 

4. Conclusion 

In this paper a uniform approximation for a broken rational torus consisting of four periodic 
orbits has been presented. The approximation was tested in a model system where the 
phase-space is mixed, and the tori are grouped around (but not too close to) a central orbit. 
It can be expected that the approximation will be even more useful in studies of globally 
near-integrable systems. 
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